Topological nematic states and non-Abelian lattice dislocations 
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An exciting new prospect in condensed matter physics is the possibility of realizing fractional 
quantum Hall (FQH) states in simple lattice models without a large external magnetic field. A 
fundamental question is whether qualitatively new states can be realized on the lattice as compared 
with ordinary fractional quantum Hall states. Here we propose new symmetry-enriched topological 
states, topological nematic states, which are a dramatic consequence of the interplay between the 
lattice translation symmetry and topological properties of these fractional Chern insulators. When 
a partially filled flat band has a Chern number N, it can be mapped to an iV-layer quantum Hall 
system. We find that lattice dislocations can act as wormholes connecting the different layers and 
effectively change the topology of the space. Lattice dislocations become defects with non-trivial 
quantum dimension, even when the FQH state being realized is by itself Abelian. Our proposal 
leads to the possibility of realizing the physics of topologically ordered states on high genus surfaces 
in the lab even though the sample has only the disk geometry. 
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I. INTRODUCTION 

Some of the most important discoveries in condensed 
matter physics are the integer and fractional quantum 
Hall (IQH and FQH) states 1-4 , which provided the first 
examples of electron fractionalization in more than one 
dimension and paved the way for our current understand- 
ing of topological order 5-7 . Conventionally realized in 
two-dimensional electron gases with a strong perpendic- 
ular magnetic field, these states exhibit a bulk energy 
gap and topologically protected chiral edge states. The 
topological order of the FQH states is characterized by 
topology-dependent ground state degeneracies 8 and frac- 
tionalized quasi-particles 4 , while the quantized Hall con- 
ductance is determined by a topological invariant - the 
Chern number - which for a band insulator (IQH states) 
can be determined by the momentum-space flux of the 
Berry's phase gauge field. 9 . 

Since Chern numbers are generic properties of a band 
structure, it is natural to expect that quantum Hall 
states can also be realized in lattice systems without 
an applied magnetic field 10,11 , and indeed such "quan- 
tum anomalous Hall" (QAH) states may be realizable 
experimentally 12 ' 13 . Recently, numerical evidence of 
fractional QAH (FQAH) states - equivalent to the 1/3 
Laughlin state and non-Abelian Pfaffian state - has 
been found in interacting lattice models with (quasi-)flat 
bands with Chern number C = l 14-21 . A flat band with 
C = 1 is similar to a Landau level, where the kinetic en- 
ergy is quenched and interaction effects arc maximized. 
Recently, a systematic wave function approach was in- 
troduced that demonstrates how to associate each FQH 
state in a Landau level to a counterpart in a generic 
C = 1 band 22 . 

While known FQH states can be realized on the lat- 
tice (see also 23-26 ), a fundamental question is whether 
new states can emerge as well. Here we show that 
the answer is yes. There are two key reasons that the 
FQAH system is different: the possibility of a band with 
C > 1, and the interplay of topological order with lattice 



symmetries 6,27,28 . While each Landau level has C = 1, a 
single band in a lattice system can in principle carry ar- 
bitrarily high Chern number. For example a model with 
C = 2 quasi-flat bands has been recently proposed 29 . 
Here we show that there are new topologically ordered 
states - topological nematic states - in a partially filled 
band with C > 1. These states carry a nontrivial repre- 
sentation of the translation symmetry and spontaneously 
break lattice rotation symmetry 30 . By using the Wannicr 
function representation 22 , a band with Chern number C 
can be mapped to C layers of Landau levels, but the C 
layers are cyclically permuted under certain lattice trans- 
lations. This leads to a dramatic interplay of these states 
with the lattice translation symmetry: a pair of lattice 
dislocations connecting different layers corresponds to a 
"wormholc" in the C-layer FQH system. The lattice dis- 
locations effectively change the topology of the space, 
giving rise to topological degeneracies. Surprisingly, the 
topological degeneracy associated with dislocations is 
nontrivial even when the state itself, in the absence of 
dislocations, is an Abelian topological state 31-34 . Our re- 
sult provides a new possibility of realizing exponentially 
large topological degeneracies without using a "genuine" 
non-Abelian FQH state, and effectively provides a way to 
experimentally observe the topology-dependent ground 
state degeneracies of FQH states. 



Wannier function description of Chern insulators 

We begin by reviewing the one-dimensional Wannicr 
state description of Chern insulators 22,35 . For a band 
insulator with N bands, the Hamiltonian is given by 



H = d 2 k4/i(k)c k 



(1) 



with /i(k) a, N x N matrix and the annihilation operator 
Ck an n-component vector. In the current work we will 
focus on the systems with only one band occupied. De- 
noting the occupied band by |k), the Berry phase gauge 
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field is cii(k) = — i (k| d/dki |k). The one-dimensional 
Wannier states are defined by 
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The phase choice ip(k x ,k y ) is determined by the con- 
dition that the states be maximally localized in the x- 
direction, which is determined by the eigenvalue equation 
x\W(k y ,n)) = x n \W(k y ,n)), where x is the x-direction 
position operator projected to the occupied band 36 . The 
formalism here can be generalized to a torus with finite 
L x ,L y , in which case the Wilson line integral is replaced 
by a discrete summation 37 . 

The Wannier states form a complete basis of the oc- 
cupied subspace, and each Wannier state is localized in 
the x direction and a plane-wave in the y direction. The 
essential property of the Wannier state is that its center 
of mass position x n (k y ) = (W(k y ,n)\x\W(k y ,n)), which 
is also the eigenvalue of the projected position operator 
x, is determined by the Wilson loop 
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Consequently, the Chern number 9 of the band C\ = 
j- J dk x dk y (d x a y — d y a x ) is equal to the winding num- 
ber of x n (k y ) when k y goes from to 27r: 



dx n {ky) 
dku 



dk v = C\ 



(5) 



In other words the Wannier states satisfy the twisted 
boundary condition \W (k y + 2tt , n)) = \W(k y , n + C\)) 
and correspondingly x n (k y + 2tt) = x n (k y ) + C\. Due to 
such a twisted boundary condition, the Wannier states 
for C\ = 1 can be labeled by one parameter K = 
k y + 27m, and the Wannier state basis \Wk) is in one- 
to-one correspondence with the Landau level wavcfunc- 
tions in the Landau gauge in the ordinary quantum Hall 
problem 22 . For C\ = 2, we have \W (k y + 2tt , n)) = 
\W(k y , n + 2)), so that the Wannier states on even and 
odd sites are two distinct families, as shown in Fig. 1 (a). 
By adiabatic continuation of the momentum k y one can 
define 
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: \W{k y ,2n-l)) 
= \W(k y ,2n)) 
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FIG. 1: A. Each state shifts over two lattice spacings: {£) — > 
(£■) +2, as k y — > ky + 2n. Thus there are two families of states 
(red and blue lines). B. The states can be mapped to two 
families of states, each parametrized by a single parameter 
K y . C. Illustration of the fact that the Chern number 2 lattice 
system is mapped to a bilayer quantum Hall system, with the 
two layers corresponding to the two families of Wannier states 
shown in panel B. 



such that 



• K J arc both continuous in the parameter 

K, and for the same K they are related by a transla- 
tion in the x direction. Compared with the C\ = 1 case, 
one can see that each family of states \ W] ( ) or |M^|-) is 
topologically equivalent to the Wannier states of a C\ = 1 
Chern insulator, which is in turn topologically equivalent 
to a Landau level quantum Hall problem. Therefore the 
Wannier state representation when C\ — 2 defines a map 
from the C\ = 2 Chern insulator to a bilayer quantum 
Hall problem, with the Wannier states on the even and 
odd sites mapped to two layers of Landau levels. In con- 
trast to a typical bilayer quantum Hall system, the two 
"layers" are now related by a lattice translation. This 
is the key observation which leads to the topological ne- 
matic states that we will propose below. 

Multilayer FQH states can exhibit a rich variety 
of different topological states; recently, a large clas- 
sification of possibilities was developed and applied 
to two-component FQH states. 38 ' 39 The simplest two- 
component generalizations of the Laughlin states are the 
(mill) states with the wavefunction 40 

$ ({Zi} , {wi}) = Y[ 0* - zj) m (wi - Wj) n Y[ (Zi - Wjf 



exp 



hi 



/4l 



(7) 



where Ib is the magnetic length, Zi = Xi + iyi are the 
compex coordinates of the ith particle in one layer, and 
similarly u>i are the complex coordinates for the other 
layer. These states can be written down for the frac- 
tional Chern insulators by simply passing to the occu- 
pation number basis, ^({nj}), where n\ is the occupa- 
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tion number of the ith orbital associated with the Ith 
layer. There are also intrinsically multilayer non-Abelian 
states 38,39 ' 41,42 . In the current paper we will focus on the 
(mml) states, where m ^ I for incompressible states. 



II. INTERPLAY WITH LATTICE 
TRANSLATION SYMMETRY AND 
DISLOCATIONS 

A 2D lattice is invariant under two independent trans- 
lation operations T x and T y . Their action on the Wannier 
states defined in Eq. (2) and (6) is 

T x \Wk) = \W%), T x \W] c ) = \Wk + 2„), 
T y \W a K ) = e iK \W a K ). (8) 

Thus T x exchanges the two sets of Wannier states but T y 
does not. 

Now consider the effect of dislocations 53 ; these are 
characterized by a Burgers vector b, which is defined 
as the shift of the atom position when a reference point 
is taken around a dislocation 43 . An x dislocation with 
b = x is illustrated in Fig. 2 A. Far away from a dis- 
location, the lattice is locally identical to one without 
a dislocation, so the dislocation is, as far as the struc- 
ture of the lattice is concerned, a point defect. Now 
consider a bilayer (mmn) state realized on the lattice 
with a dislocation. As is shown in Eq. (8), the two 
sets of Wannier states are related by translation in x di- 
rection. Thus when one goes around an x-dislocation, 
the two "layers" consisting of Wannier states | ) and 
|Wjf) are exchanged. The map defined by the Wannier 
states, which maps the C\ = 2 Chern insulator to a bi- 
layer FQH system, maps the Chern insulator on a lattice 
with a pair of dislocations to a bilayer FQH state defined 
on a "Ricmann surface" with a pair of branch-cuts, as 
is illustrated in Fig. 2 B. This is the key observation 
which indicates that the x-dislocations in this system 
have nontrivial topological properties. By comparison, 
the y-dislocations do not exchange the two layers and 
thus do not correspond to a topology change in the ef- 
fective bilayer description. 

III. TOPOLOGICAL DEGENERACY OF 
DISLOCATIONS 

Although the (mml) quantum Hall state considered is 
Abclian, the ^-dislocation carries a nontrivial topological 
degeneracy. 54 To understand this, start from the simplest 
case of (mmO) state, which is a direct product of two 
Laughlin states. For such a state, the Chern insulator on 
a torus is mapped to two decoupled tori with a Laugh- 
lin 1 jm state defined on each of them, with total ground 
state degeneracy of m? . When a pair of x-dislocations 
is introduced, the two tori arc connected by the branch- 
cut. If wc do a reflection of the top layer according to 




FIG. 2: A: Illustration of an x-dislocation. B: (Upper pan- 
nel) Illustration that an x-dislocation leads to a branch cut 
around which the two effective layers are exchanged. (Lower 
pannel) A reflection of the top layer maps the branch cut be- 
tween a pair of dislocations into a "worm hole" connecting 
the two layers. C: A torus with two pairs of s-dislocations is 
equivalent to two tori connected by two "worm holes" , which 
is a genus 3 surface. This picture illustrates the fact that 
dislocations carry nontrivial topological degeneracy. 



the x axis, the branch cut becomes a "worm hole" be- 
tween the two layers, as is illustrated in Fig. 2 B. Thus 
the two tori are connected, resulting in a genus 2 sur- 
face. For two pairs of dislocations, the two layers are 
connected by two worm holes and the whole system is 
topologically equivalent to a single Laughlin 1/m state 
on a genus 3 surface, as is shown in Fig. 2 C. Thus the 
ground state degeneracy becomes m 3 . In general, when 
there are 2n rr-dislocations on the lattice, the space is 
effectively a genus n + 1 surface and the ground state de- 
generacy for n > is m n+1 . It follows that the average 
degree of freedom carried by each dislocation-known as 
the quantum dimension-is d = y/m. Thus we can see 
that the x-dislocation carries a nontrivial topological de- 
generacy, in the same way as a non-Abelian topological 
quasiparticle. 

The discussion above can be generalized: for the (mml) 
state, n > pairs of dislocations on a torus leads to the 
topological degeneracy of \m 2 — l 2 \\m — so that the 

quantum dimension of each dislocation is d = y[m — l\ 
(recall m ^ I for incompressible FQH states). For I ^ 0, 
the system cannot be mapped to a Laughlin state on 
high genus surface as the two layers are not decoupled 31 . 
The topological degeneracy can be computed from the 
bulk Chern-Simons effective theory 31 . In the following we 
provide an alternative understanding of the topological 
degeneracy using the edge states, as it is more rigorous 
for I =/= 0, and it helps to provide a clearer understanding 
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of the topological degeneracy. 



IV. TOPOLOGICAL DEGENERACY FROM 
EDGE STATE PICTURE 

Here we will study in detail the topological degeneracy 
in the topological nematic states with dislocations by us- 
ing an edge state picture. By cutting the FQH state on a 
compact manifold along a line, one obtains a FQH state 
with open boundaries and gapless counter-propagating 
chiral edge states on the boundary. The FQH state be- 
fore the cut can be obtained by "gluing" the edge states 
back by introducing intcr-edge electron tunneling (Fig. 
3A). The topological degeneracy comes from the fact 
that there are in general multiple degenerate minima 
when the electron tunneling becomes relevant 44 . In other 
words, the topologically degenerate ground states of the 
bulk topological system correspond to degenerate ground 
states due to spontaneous symmetry breaking in the edge 
theory. 

The topological degeneracy of the x-dislocations is in- 
dependent of their location, so to compute the topologi- 
cal degeneracy we may orient the n pairs of x-dislocations 
along a single line. Next, we cut the system along this 
line to obtain a gapless CFT along the cut, and then 
we understand the resulting ground state degeneracy by 
coupling the two sides of the cut in the appropriate way. 

Below, we first review the chiral edge theory of Abelian 
FQH states and how to understand their torus degener- 
acy from the point of view of the edge theory. Next, 
we show how this cut-and-glue procedure can be used to 
compute the ground state degeneracy in the presence of 
dislocations for general two-component Abelian states. 

The edge theory for an Abelian FQH state described 
by a generic A"-matrix is given by the action 44,45 



dye 



1 
4?T 



dxdt[Kud t (l)Lidx(f>Lj - Vud x (f> L id x (j) 



LJ 



(9) 



where Hi denotes left-moving chiral bosons for I = 
1, • • • , dim K . Here and below the repeated indices 7, J 
are summed. The field Hi is a compact boson field with 
radius R = 1: 

Hi ~ Hi + 27r - (10) 

Quantizing the theory in momentum space yields 6 

[d x Hi(x), Hj{v)} = 2nKj J 1 6(x - y). (11) 

Integrating the above equation gives: 

[<j>Ll( x )i <pL.j(y)} = TrKyJsgn(x - y) (12) 

The electric charge density associated with Hi is given 
by 



Pli = ^-d x (f)Li> 
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FIG. 3: The edge state understanding of the topological de- 
generacy. A. The dislocations are oriented along a single line, 
and then the system is cut along the line, yielding gapless 
counterpropagating edge states along the line. The original 
FQH state is obtained from gluing the the system back to- 
gether by turning on appropriate inter-edge tunneling terms. 
B. Depiction of the two branches (red and blue) of counter- 
propagating edge excitations. The arrows between the edge 
states indicate the kinds of electron tunneling terms that are 
added. Away from the dislocations, in the A regions, the 
usual electron tunneling terms involving tunneling between 
the same layers, ^\ RI ^ e Li + H.c, are added. In the regions 
including the branch cuts separating the dislocations, twisted 
tunneling terms are added: ^f eR1 ^ e L2 + ^ e m ^eli + H.c. oti 
indicate the mid-points of the A or B regions. 



and the 7th electron operator is described by the vertex 
operator 



^eLi = e lKlj4 > LJ . 



(14) 



Note that normal ordering will be left implicit (ie 
e iKjj<t> LJ EE . e iKu4>Lj jf we consider the FQH state 
on a cylinder, we will have a left-moving chiral theory on 
one edge, and a right-moving chiral theory on the other 
edge. For the right-moving theory, the edge action is 

Sedge = ^ / dxdt[-K ud t (j)Rld x (j) R J - Vud x (f> R A<t>Rj}, 

(15) 

the charge is 



Pri = —d x <j)Ri, 

Z7T 



and the electron operator is 



ri = e 



-iKr 



(16) 



(17) 



(13) 



Now if we bring the edges close together, the electrons 
can tunnel from one edge to the other. The electron 
tunneling operators are: 

lY.^eLI^eRI + *t RI *eLl] = ^ h COS^jj^), 

(18) 
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where 

<i>J = <t>LJ + 4>RJ 

is a non-chiral boson. Note that 

[fa( x ),fa( y )] = o. 



(19) 



(20) 



The tunneling terms gap out the edge and lead to a set of 
degenerate minima. Assuming that tj < 0, the minima 
occur when 



K 



U9J 



(21) 



where pi is an integer. That is, the minima occur when 
fa- = 2-kKJjPj. Note that since fa- ^ fa- + 2n, the 
degenerate states can be labelled by an integer vector p, 
which denotes the eigenvalues of of e 1 ^ 1 : 



(p\e»'\p) 



(22) 



for integers pj . Two different vectors p and p' arc equiv- 
alent if taking (f>j — > fa- + 2Trnj, for integers nj takes 
p —> p + Kp = p' . Thus the ground state degeneracy 
on a torus is given by the number of incquivalcnt integer 
vectors p, which is determined by |Det K\. 46 

For what follows, we specialize to the case where K 
m I 



is a 2 x 2 matrix, K 



I 



On a torus, there are 



|Det K\ = \(m — l)(m + l)\ different ground states. From 
the edge theory, these can be understood as cigenstates 
of e *(4>i±<fe) ; wit h eigenvalue e 27Tt P±^ m±l \ where p± are 
integers. Therefore, the eigenstates can be labelled by 
(P+,P-), where 

((p + ,p_)|e l ^ 1± ^)|(p + ,p_)) =e 2**P±/(m±0 (23) 

Now suppose we have n pairs of dislocations, such that 
going around each dislocation exchanges the two layers, 
as discussed in the text. Each pair of dislocations is sep- 
arated by a branch cut. Let us align all the dislocations, 
and denote the regions without a branch cut as Ai, and 
the regions with a branch cut as Bi (Fig. 3A). Now imag- 
ine cutting the system along this line, introducing coun- 
terpropagating chiral edge states. The gapped system 
with the dislocations can be understood by introducing 
different electron tunneling terms in the A and B regions 
(Fig. 3B): 



) ^\ L1 ^ eR i + *L 2 * efl 2 + H.c if X G M 



2 I *eLl*ei?2 + ^[ L2 ^ ' eRl + H.C if X £ B t 

(24) 



Introducing the variables 



5L1 

h L 2 



»R2 
t>Rl 



we rewrite (24) as 



Y1i<x>b( k ij<Pj) if x e Ai 
^2jCos(K u4>j) if x G Bi 



(25) 



(26) 



We can write the ground state approximately as a ten- 
sor product over degrees of freedom in the different re- 
gions: 



(27) 



where |a.;) is a state in the Hilbert space associated with 
Ai, and similarly for Bi. The set of ground states is 
a subspacc of the |Det K\ 2n states associated with the 
minima of 8C tun n for each region. 

If the interaction between different regions are ignored, 
one would have obtained \m 2 — l 2 \ degenerate states for 
each region labeled by cigenstates of (<j>i(x) ± (j>2(x)) or 
(4>\{x) ±4>2(x)) in each region. However, we observe that 
phases in A and B regions cannot be simultaneously di- 
agonalizcd. Instead, 



<j>i{x),(j>j(x') 



-1 1 ) sgn ^ ~ x "> 
J ij 

(28) 



Therefore 



e ir t>i(x) e «Oi<y)-02O')) e -i<t>i(x) — e i ^5l s g n ( x - x ') e i (0i( a: ')-'/ , 2(a;')) 



J4>i(x) e i(4>i(x')+<f> 2 (x')) — Q 



(29) 



From these commutation relations we see that </>i + 4>2 
can have common eigenstates with tfii and <j>2 in the A 
regions, but the eigenstates of 4>i and fa are necessar- 
ily a superposition of the m — I eigenstates of fa — fa. 
From this consideration, we see that we can pick a basis 
where |Det K\ n possible states are associated with the A 
regions, and the B regions contribute (m + l) n possible 
states, for a total 

iV = |Dct K\ n (m + l) n = (m 2 - l 2 ) n {m + l) n (30) 

possible states. 

However, there are also 2n — 1 global constraints that 
must be satisfied. The charge at each dislocation should 
be a definite quantity, which, without loss of generality, 
we set to 0. This is because states at the dislocation 
which carry different charge cannot be topologically de- 
generate, since charge is a local observable. These con- 
straints are expressed as 



d x {fa + fa)dx = 0,i = l,2,..,2n (31) 



where a.2%-1 and ct^i are the mid points of region Ai and 
G Bi, respectively, so that the integration region includes 
the «-th dislocation (Fig. 3B). By definition £\ Qi = so 
that the number of independent constraints (for n > 0) is 
2n— 1. Each constraint reduces the number of states by a 
factor of m + 1 since only the sector of fa + fa is involved. 
Therefore the topological degeneracy in the presence of 
n > pairs of dislocations is 



(m + l) 2n ~ 



= Urn 2 - l 2 )(m - l) r 



(32) 
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V. MORE GENERIC TOPOLOGICAL 
NEMATIC STATES AND EFFECTIVE 
TOPOLOGICAL FIELD THEORY 

As shown above, the topological degeneracy is only as- 
sociated with the ^-dislocations and not to ^-dislocations. 
Such states apparently break rotation symmetry but pre- 
serve translation symmetry, which is why we name them 
topological nematic states. The reason for the disparity is 
that the FQH states are constructed using Wannicr states 



W 



K 



localized in the ^-direction. If we instead define 



the mapping from the lattice system to the bilaycr FQH 
system using Wannier states localized in the y-direction, 
the resulting state will have topological degeneracy asso- 
ciated with ^-dislocations. These two different types of 
topological states can be related by redefining the Bril- 
louin zone. More generally, we can choose the reciprocal 
lattice vectors to be ei = (p,q) and e,2 = (p',q') with 
|ex x e2| = pq' — qp' = 1, and define one-dimensional 
Wannier states by the Fourier transform of Bloch states 
along the ei direction, as shown in Fig. 4. The Wannier 
states obtained in this way are localized along ei direc- 
tion. For a dislocation with Burgers vector b = (b x ,b y ), 
around a dislocation the translation along the ei direc- 
tion is b ■ ei = b x p + b y q. Therefore the exchange of 
two layers and topological degeneracy only occurs when 
b ei = b x p+b y q = 1 mod 2. This condition classifies the 
topological nematic states through four types defined by 
(p, q) mod 2 = (0, 0), (0, 1), (1, 0), (1, 1). The topological 
degeneracy is associated with x (y)-dislocations if and 
only if p (q) is odd. Rigorously speaking, (0, 0) cannot 
be realized by any state constructed in this way, because 
pq' — qp' = 1 requires that at least one of p and q is 
odd. However one can interpret the type (0, 0) as the 
ordinary bilayer FQH state, where the two layers are not 
exchanged under any translation. 

To describe the topological nematic states better, we 
have developed an effective field theory that naturally 
incorporates the interplay between topological proper- 
ties and the lattice dislocations in this system. Note 
that a dislocation in a two-dimensional crystal is a defect 
with long range interaction, similar to superfluid vortices. 
Thus the effective field theory is different from the topo- 
logical field theory studied in Ref. 31 which describes the 
phase with finite energy twist defects. The proper effec- 
tive field theory should satisfy the following conditions: 

1. The theory describes the dynamics of the crystal, 
which is characterized by the displacement field 
u = {u Xl Uy). In a 2d crystal the translation sym- 
metry Rx R is broken to Z x Z, so that the order 

parameter u is periodic + 1, Uy ~ Uy + 1. 

Therefore u G U(l) x U(l). 

2. The theory should reduce to a U(l) x U(l) Chern- 
Simons theory in the absence of dislocations. 

3. When dislocations are present, the two U{1) gauge 




FIG. 4: Illustration of different definitions of Wannier states, 
which lead to different types of topological nematic states, ei 
and e2 are two reciprocal vectors defining a Brillouin zone. 
The Wannier state basis can be constructed by taking the 
Fourier transform of Bloch states along one periodic direction 
of the Brillouin zone, marked by the red, blue and purple lines 
with arrows. The red, blue and purple lines correspond to 
topological nematic states of the type (1, 0), (0, 1) and (1, 1), 
respectively (see text). 



groups should be exchanged when we take a refer- 
ence point around a dislocation point. 

A natural effective theory satisfying the conditions 
above can be constructed by using a U(2) gauge field 
dp coupled to a Higgs field H of the form H = a ■ ne t8 ( u \ 
with a are the three Pauli matrices, n a real unit vector, 
and 



6(u) 



u ■ ei 



(33) 



is a phase factor determined by the displacement field u 
. Here ei is the reciprocal vector defining the Wannier 
states as is illustrated in Fig. 4 and earlier in this sec- 
tion. H is a traceless, unitary 2x2 matrix in the adjoint 
representation of t/(2), which transforms under a gauge 
transformation H —> g~ 1 Hg, with j a (7(2) matrix, n 
transforms as a vector under the SU (2) subgroup of the 
U(2) and e 10 is gauge- invariant. The effective Lagrangian 
has the following form: 



C 



l r 



p(d t u) - XijkidiUjd k ui 



I 



-in 



-e^ T Tr 



2i 

O 



4?r 



e^ T Tr [o„] d v Tr [a T ] + JTr [D^D^H} (34) 



where the covariant derivative is D^H = d^H + i [a 



The first term is the standard elasticity theory of the 
crystal. The Chern-Simons fields describe the topological 
degrees of freedom of the electrons, which is coupled with 
the crystal through the Higgs field H . 

We first consider a system without dislocation. For 
example for u — one can choose a constant Higgs 



field H = <r z , which gives mass to two of the four 
components of the U{2) gauge field a. The part of 
that remains massless is: a M = 0^5 + Ou%-- Denoting 



u(d) 



h (a° ± al) , the CS term for 



and 



reduces 



to 



C CS = T ^ VT (mald v a u T + mafaaf + 2la^a d T ) 



(35) 



which correctly recovers the U (1) x U(l) CS theory of the 
{rami) state 6 ' 47 , with u and d labeling the two layers. 

Now we consider the effect of dislocations. Around a 
dislocation with Burgers vector b, changes by 7rei • b. 
Therefore if ei • b = 1 mod 2, e lB has a half wind- 
ing around the dislocation. This is allowed if and only 
if n also has a half-winding, compensating for the mi- 
nus sign from the shift of 6. Such a topological de- 
fect is similar to the half vortex in a spinful {p + ip) 
superconductor. 48 Different from the p + ip superconduc- 
tor in which the SU(2) charge is global and the U(l) 
charge is coupled to the electromagnetic gauge field, in 
the current case the SU (2) charge carried by the vector 
n is gauged and the U{1) charge remains global. Math- 
ematically, the space of the order parameter described 
by H is S 2 x U{l)/Z 2 ; the Z 2 corresponds to identifying 
e lS n with e l ^ 6+ ^ {— n). This has a nontrivial fundamen- 
tal group: 7Ti (S 2 x U{l)/Z 2 ) = Z 2 , allowing topologi- 
cal^ non-trivial point defects. The invariant subgroup of 
U (2) preserving the Higgs field is U{1) xU(l), defined by 
the rotations g = e l ( a +^ n CT )/ 2 . When we adiabatically 
take a point around a dislocation, the vector n(r) is adi- 
abatically rotated to — n(r), such that the two U{l)'s in 
the invariant subgroup are exchanged. Thus the effective 
theory (34) reproduces the fact that the two U{1) gauge 
fields, describing the two layers, in the U{l)xU{l) Chern- 
Simons theory are exchanged when a reference point is 
taken around a dislocation. By construction, this theory 
also correctly describes the fact that only dislocations 
with ei • b odd have topological degeneracy. Due to the 
gapless Goldstonc modes described by displacement field 
u, the dislocations have logarithmic interaction. 



VI. EVEN/ODD EFFECT AND DETECTION IN 
NUMERICS 



Even in the absence of dislocations, the topological 
nematic states exhibit unique topological characteristics 
that can aid in their detection in numerical simulations. 
To see this, below we will focus on the {rami) topological 
nematic states of type (1,0), defined on a lattice with 
periodic boundary conditions. 

Let us first consider the situation shown in Fig. 5 (a), 
where we suppose that along an entire line parallel to 
the x-direction, the hoppings have been twisted. We also 
suppose that the number of sites in the x-direction is 
even. Conceptually, we can imagine constructing this by 
taking a pair of dislocations, moving one of them around 




B 



2n + 1 



FIG. 5: Illustration of two lattice configurations to detect 
the topological nematic states, (a) A torus with a kink in 
hopping along the dash line, and with even number of sites 
in x direction, (b) A torus with no kink in hopping, but with 
odd number of sites in x direction. In both cases, periodic 
boundary condition is imposed to both directions. For both 
configurations, the {mml) topological nematic state of type 
(1, 0) has a reduced ground state degeneracy of |m + Z| instead 
of the degeneracy of |m 2 — 1 2 \ on a torus with even number of 
sites in x direction. 



the x-direction of the torus, and re-annihilating them. As 
we explain below, this change in the geometry effectively 
changes the topology, as in the case with dislocations. To 
see this, consider the (mmO) states, which can be viewed 
as a direct product of two Laughlin 1/m states. For a 
regular lattice with periodic boundary conditions, the low 
energy theory can be mapped onto two decoupled tori, 
each giving a degeneracy of m, for a total topological 
ground state degeneracy of m 2 . However, the twisted 
hopping has the effect of gluing the two tori together 
into a single torus, and so we expect that the ground 
state degeneracy is actually m. We can further see this 
in an alternative way by using the edge theory picture as 
presented earlier. By introducing twisted couplings into 
the edge theory, we introduce the tunneling terms 



6S, 



twisted 



^2 C0S ( K uh), 



(36) 



with K 



I 



I 



and 4>j defined in Eq. (25). Naively, 



this would give |DctAT| different states labeled by the 
eigenstates of e 1 ^ 1 ^ = e lp± ^ m±l \ similar to the dis- 
cussion in Sec. IV. However, we observe that there is a 
gauge symmetry associated with the independent fermion 
number conservation within each layer, which is associ- 
ated with the shift symmetry 



VLi 



T>Li 



(37) 



where 6\ and 6 2 are independent. In the absence of tun- 
neling between the different layers, this fermion number 
conservation is clearly present. Since the {mml) state is 
gapped, tunneling between the different layers is an ir- 
relevant perturbation. Therefore even in the presence of 
tunneling between the different layers, we may regard the 
fermion number conservation symmetry within each layer 
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as an emergent symmetry of the topological state. We 
conclude that, of the \(m + l)(m — l)\ degenerate minima 
of SStwisted, the ones associated with the same eigenvalue 
f e «0i+202 an( f different eigenvalues of e z ^ 1_i< ^ 2 arc gauge 
equivalent to each other. Therefore, we conclude that the 
topological degeneracy of such a twisted lattice is |m + /|. 
In contrast, in the case of the non- twisted hoppings the 
gauge symmetry acts in a trivial way on all the states 
and therefore does not affect the topological degeneracy 
|Detif|. 

Alternatively, consider a regular lattice, but with an 
odd number of sites along the x-direction, as shown in 
Fig. 5 (b). In this case, there will be a line at some point 
in x where Wannicr functions associated with different 
layers are adjacent. Cutting the system along this ver- 
tical line, we obtain gapless edge states, and we observe 
that in gluing the system back together, we again have 
twisted hoppings along the entire vertical line. Following 
the same argument as the last paragraph, in this case 
the ground state degeneracy is |m + 1\. This leads to a 
remarkable signature of the topological nematic states. 
For a regular lattice, when the number of sites in the x- 
direction is even, the topological degeneracy is |m 2 — l 2 \. 
When it is odd, the topological degeneracy is \m + l\. 
These topological degeneracies are independent of the 
length in the y-dircction. Such an even-odd effect can 
easily be used as a numerical diagnostic to test for such 
topological nematic states. 



O T t 



FIG. 6: (a) Illustration of the domain wall between two topo- 
logical nematic states (1,0) (upper half) and (0,1) (lower 
half), (b) Illustration that the left and right moving edge 
states along the domain wall transform differently in lattice 
translation T x . The right movers from the (1,0) state are ex- 
changed by T x while the left movers from the (0, 1) state are 
separately translation invariant. 



tunneling terms must involve the symmetric combination 
cos (4>li + 4>L2 + 4>Ri + 4>R2)- Therefore, either transla- 
tion symmetry is preserved and the mode associated with 
the relative combination 4>\l + 4>ir — 4>2L — 4>2R is gapless, 
or translation symmetry is spontaneously broken along 
the domain wall. This provides a novel example of trans- 
lation symmetry protected gapless edge states. 



VIII. DISCUSSION AND CONCLUSIONS 



VII. DOMAIN WALLS AND TRANSLATION 
SYMMETRY-PROTECTED GAPLESS MODES 

If the strongly interacting Hamiltonian that realizes 
the topological nematic states has a discrete lattice rota- 
tion symmetry, then the topological nematic states asso- 
ciated with the class (1,0) or (0,1) (see Section V) spon- 
taneously break this discrete rotation symmetry. In such 
a situation, a physically realistic system can be expected 
to consist of domains involving different orientations. At 
the domain walls, it is possible that translation symmetry 
is preserved along the domain wall. As we show below, 
the translation symmetry along the domain wall can pro- 
tect a single bosonic channel from being gapped out. 

For concreteness, consider a domain wall between a 
(0, 1) and a (1,0) topological nematic state, as shown in 
Fig. 6. Before coupling the two domains together, each 
one has two branches of chiral gapless edge states: <f>Li 
and 4>m, for i = 1, 2, where L/R stand for the left /right- 
movers. Naively, these counterpropagating modes will 
be gapped out by inter-edge tunneling. However, they 
have different properties under translation symmetry. In 
particular, under a translation along the domain wall, 

4>L1 -> <f>LU <t>L2 ™> 

<t>Rl -> 4>R2, 4>R2 -> <t>Rl, (38) 

We see that for the translation symmetry along the do- 
main wall to be preserved, the only allowed inter-edge 



This realization that lattice dislocations carry non- 
trivial topological degeneracy raises a host of possible 
new directions. Conceptually, we are now in need of 
a more general theory of topological degeneracy of dis- 
locations and their possible braiding properties in gen- 
eral FQH states, both in bands with higher Chcrn num- 
bers and in non-Abelian FQH states. For example, for a 
Chern insulator with Chern number C\ = N, the same 
construction maps the lattice system to a iV-layer FQH 
state. Dislocations are "branch cuts" with degree N 
in such a system, around which the N layers are cycli- 
cally permuted. It would also be important to verify 
our predictions through numerical studies. The even- 
odd effect discussed in Sec. (VI) is expected to generalize 
to a mod N behavior in the topological degeneracy for 
d = N. 

Note that while the dislocations studied here have non- 
Abelian properties, they are confined excitations of the 
system because the energy cost of separating dislocations 
grows logarithmically with their distance. Nevertheless, 
physical materials can easily have many lattice disloca- 
tions, and if FQH states arc realized in a material with 
higher Chern bands, we expect that the lattice disloca- 
tions, in the limit of low dislocation density, can induce 
an extensive topologically stable entropy for tempera- 
tures T > Ae^ 1 ^, where £ cx 1/A is the correlation 
length of the gapped FQH state, A is the energy gap, 
and / is the typical spacing between dislocations. 

A further interesting, exotic possibility is to imagine 
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quantum melting the lattice to deconfine these disloca- 
tions. The resulting states, if realizable, are "topological 
liquid crystals" and may correspond to the non-Abelian 
orbifold states proposed in 49 and that are described by 
U(l) x U{1) x Z 2 CS theory. 31 

Acknowledgements 

We would like to acknowledge inspiring discussions 
with Alcxci Kitaev and Xiao-Gang Wen and thank the 



Simons Center for Geometry and Physics, Stony Brook 
University for hospitality where these discussions were 
made. We would like to acknowledge helpful discussions 
with Zheng-Cheng Gu, Steven A. Kivelson and Cenke 
Xu. This work was supported by a Simons Fellowship 
(MB) and the Alfred P. Sloan foundation and the David 
and Lucile Packard Foundation (XLQ). We thank the 
KITP program Topological Insulators and Superconduc- 
tors for hospitality while part of this work was completed. 
This research was supported in part by the National Sci- 
ence Foundation under Grant No. NSF PHY11-25915. 



1 K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 
45, 494 (1980). 

2 D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. 
Lett. 48, 1559 (1982). 

3 R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 

4 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

5 S. C. Zhang, Int. J. Mod. Phys. B 6, 25 (1992). 

6 X.-G. Wen, Quantum Field Theory of Many- Body Systems 
- From the Origin of Sound to an Origin of Light and 
Electrons (Oxford Univ. Press, Oxford, 2004). 

7 C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. D. 
Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

8 X. G. Wen, Phys. Rev. B 41, 12838 (1990). 

9 D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 
M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

10 F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

11 X. L. Qi, Y. S. Wu, and S. C. Zhang, Phys. Rev. B 74, 
085308 (2006). 

12 C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, 
Phys. Rev. Lett. 101, 146802 (2008). 

13 R. Yu, W. Zhang, H. J. Zhang, S. C. Zhang, X. Dai, and 
Z. Fang, Science 329, 61 (2010). 

14 K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, Phys. Rev. 
Lett. 106, 236803 (2011). 

15 N. Regnault and B. A. Bernevig (2011), arXiv:1105.4867. 

16 Y.-L. Wu, B. Bernevig, and N. Regnault (2011), 
arXiv:1111.1172. 

17 Y.-F. Wang, H. Yao, Z.-C. Gu, C.-D. Gong, and D. Sheng 
(2011), arXiv:1110.4980. 

18 T. Neupert, L. Santos, C. Chamon, and C. Mudry, Phys. 
Rev. Lett. 106, 236804 (2011). 

19 E. Tang, J.-W. Mei, and X.-G. Wen, Phys. Rev. Lett. 106, 
236802 (2011). 

20 D. N. Sheng, Z.-C. Gu, K. Sun, and L. Sheng, Nat. Comm. 
2, 389 (2011). 

21 D. Xiao, W. Zhu, Y. Ran, N. Nagaosa, and S. Okamoto, 
Interface engineering of quantum hall effects in digital het- 
erostructures of transition-metal oxides, arXiv:1106.4296 
(2011). 

22 X.-L. Qi, Phys. Rev. Lett. 107, 126803 (2011). 

23 J. McGreevy, B. Swingle, and K.-A. Tran (2011), 
arXiv:1109.1569. 

24 S. Prameswaran, R. Roy, and S. Sondhi (2011), 
arXiv:1106.4025. 

25 A. Vaezi (2011), arXiv:1105.0406. 

26 G. Murthy and R. Shankar (2011), arXiv:1108.5501. 

27 X.-G. Wen, Phys. Rev. B 65, 165113 (2002). 



28 Y.-M. Lu and Y. Ran (2011), arXiv:1109.0226. 

29 F. Wang and Y. Ran, Nearly flat band with chern number 
c=2 on the dice lattice, arXiv:1109.3435 (2011). 

30 S. Kivelson, E. Fradkin, and V. Emery, Nature 393, 550 
(1998). 

31 M. Barkeshli and X.-G. Wen, Phys. Rev. B 81, 045323 
(2010). 

32 H. Bombin, Phys. Rev. Lett. 105, 030403 (2010), 
URL http : //link . aps . org/doi/10 . 1 103/PhysRevLett . 
105.030403. 

33 A. Kitaev and L. Kong (2011), arXiv:1104.5047. 

34 A. Kitaev, Annals of Physics 321, 2 (2006). 

35 S. Coh and D. Vanderbilt, Phys. Rev. Lett. 102, 107603 

(2009) . 

36 S. Kivelson, Phys. Rev. B 26, 4269 (1982). 

37 R. Yu, X.-L. Qi, A. Bernevig, Z. Fang, and X. Dai, e-print 
arXiv:1101.2011 (2011). 

38 M. Barkeshli and X.-G. Wen, Phys. Rev. B 82, 245301 

(2010) . 

39 M. Barkeshli and X.-G. Wen, Phys. Rev. B 82, 233301 
(2010). 

40 B. Halperin, Helvetica Physica Acta 56, 75 (1983). 

41 E. Ardonne, F. van Lankvelt, A. Ludwig, and 
K. Schoutens, Physical Review B 65 (2002). 

42 E. Ardonne, N. Read, E. Rezayi, and K. Schoutens, Nu- 
clear Physics B 607, 549 (2001). 

43 P. Chaikin and T. Lubensky, Principles of condensed mat- 
ter physics (Cambridge University Press, 1995). 

44 X. G. Wen, Phys. Rev. B 41, 12838 (1990), URL http: 
//link . aps . org/doi/10 . 1103/PhysRevB . 41 . 12838. 

45 X. Wen, Int. J. Mod. Phys. B6, 1711 (1992). 

46 X. G. Wen, Phys. Rev. B 40, 7387 (1989), URL http: 
//link . aps . org/doi/10 . 1103/PhysRevB . 40 . 7387. 

47 X.-G. Wen and A. Zee, Phys. Rev. B 46, 2290 (1992). 

48 S. B. Chung, H. Bluhm, and E.-A. Kim, Phys. Rev. Lett. 
99, 197002 (2007). 

49 M. Barkeshli and X.-G. Wen, Phys. Rev. B 84, 
115121 (2011), URL http://link.aps.org/doi/10. 1103/ 
PhysRevB. 84. 115121. 

50 V. Juricic, A. Mesaros, R.-J. Slager, and J. Zaanen, 
Universal probes of two-dimensional topological insulators: 
Dislocation and pi-flux, arXiv:1108.3337 (2011). 

Y. Ran, Y. Zhang, and A. Vishwanath, Nature Phys. 5, 
298 (2009). 

52 Compared to the corresponding expression in Ref. 22 , we 
have added the first two terms in Eq. (3) which only de- 
pends on k y and fixes the gauge choice between the Wan- 



10 



nier states of different k y . 

Several other topological effects induced by lattice disloca- 
tions have been studied in other topological states. 50 ' 51 
A related situation has been studied in two-component 
Abelian quantum Hall states 31 where Z2 twist defects that 
have the effect of exchanging QH layers were found to be 
non-abelian quasi-particles. Subsequently, Z2 twist defects 



arising from dislocations in the toric code model was stud- 
ied in 32 , where it was also found that they are non- Abelian 
anyons. The current work is the first proposal where the 
Z2 twist defect can be realized by a real lattice dislocation 
in a realistic electronic system. 



